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Abstract 

We suggest a new approach for the automatic and fully numerical evaluation of one- 
loop scattering amplitudes in perturbative quantum field theory. We use suitably 
formulated dispersion relations to perform the calculation as a convolution of tree- 
level amplitudes. This allows to take advantage of the iterative numerical algorithms 
for the evaluation of leading order matrix elements. 
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1 Introduction 



The LHC will provide unprecedented experimental and theoretical challanges. 
Among others, the presence of many final state jets for any interesting observ- 
able, makes absolutely crucial a better control of theoretical predictions and in 
particular of higher order contributions. At present several Monte Carlo event 
generators based on exact Leading Order (LO) matrix elements are available 
such as Alpgen [lj, Helac [2j, MadEvent [3] and Sherpa [4J. These codes are 
of extreme importance for the modeling of multi-jet final states. Nevertheless, 
the main problem connected with LO calculations is the strong factorisa- 
tion/renormalisation scale dependence. In this respect it would be useful to 
have also Next-to-Leading Order (NLO) calculations for such multi-jet final 
states in order to better control the overall normalization of the theoretical 
predictions. Given the huge complexity of the calculations, many theoretical 
efforts have been recently devoted to the development of new approaches to 
the problem. 

The real part of the NLO corrections of any n-body final state, given by 
(n + l)-body tree-level matrix elements, can be computed efficiently thanks to 
automated helicity amplitudes [5] or fully numerical algorithms [6]TT] . which 
allowed the development of the above mentioned LO matrix element event 
generators. On the other hand the NLO virtual corrections for a generic 2 — > n 
process are known for n = 2 and only for some 2 — ^ 3 processes. Among the 
latter it is worth mentioning the QCD corrections to pp — > ZZZ, WWZ, 
HHH [8], pp —>■ Hjj [9] and pp —>■ tij [10]. For 2^4 processes, only the 
QCD corrections to the weak boson fusion pp — > WWjj [llj, WZjj [12] 
and the complete electroweak 0(a) corrections to e + e~ — > 4 fermions [13] 
have been calculated with standard diagrammatic approaches, supplemented 
with new reduction techniques of scalar and tensor integrals [H] . No complete 
calculation of NLO QCD corrections to pp — > 4 partons exists at present. 
Only the calculation of six-gluon scattering amplitudes has been completed 
successfully [15], even if the needed CPU time renders the calculation unsuited 
for event generation. Within QED, the calculation of n— photon amplitudes 
have been performed by several groups, with n up to six [16J. 




In order to overcome the difficulties in the calculation of virtual corrections 
to multi-leg processes, several groups have developed new (semi)numerical 
and analytical techniques. It is worth mentioning the work of Ref . [17] , based 
on the method developed in Ref. [18J, which has been used to build up the 
complete two- loop renormalization program within the Standard Model [?]. 
Pure numerical algorithms have been proposed in Ref. [20], with applications 
to e + e~ annihilations. 

The analytical methods exploit the unitarity cut method and on-shell recursion 
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relations to calculate the so-called cut-constructible part of loop amplitudes 
(see Ref. [2TJ for an exhaustive review and references). The method is partic- 
ularly suited for loop calculations in gauge theories, where the final analytical 
expressions are quite simple. An efficient algebraic technique, suitable for nu- 
merical implementation, has been suggested in Ref. [22]. The main problem of 
the unitarity cut approach is related to the so-called rational part of the am- 
plitude, for which other methods are required [23] . Very recently, in Ref. [23] 
it has been proposed an algorithm based on .D-dimensional unitarity, which 
allows the calculation of both cut-constructible and rational parts of one-loop 
scattering amplitudes. Moreover the method reduces the algorithmic complex- 
ity at the polynomial level and it is suited for numerical implementation. 

The common feature of the above mentioned techniques is the use of dimen- 
sional regularization and the decomposition of the amplitude in linear combi- 
nations of basic scalar loop integrals. 

In the present note we study an alternative way, of purely numerical nature, 
for the calculation of loop-amplitudes, exploiting the fact that we have already 
at hand efficient and well developed numerical tools for the calculation of tree- 
level scattering amplitudes. We follow the approach suggested by Veltman in 
Ref. [25j, according to which any one- loop amplitude can be written as a 
convolution of LO amplitudes by means of appropriate dispersion relations. 
The method has been used in Refs. [2"B*f2"6 r ] at the formal level to discuss 
causality for individual diagrams. With a particular way of balancing energy 
conservation in cut diagrams by adding an 'artificial particle', a Feynman 
graph can be computed by the convolution of LO amplitudes, cutting any 
two internal lines. The method can be reformulated in terms of amplitudes, 
rather than in terms of individual Feynman graphs, and in this way we can 
retain the main advantage of LO approach, namely the power law growth of 
computational complexity rather than factorial growth. 

The numerical nature of the procedure is an added bonus for the automation 
of the procedure. In the calculation of a scattering cross sections, the loop 
integration is performed through the Monte Carlo method together with the 
phase space integration. This leaves two important open problems: the sta- 
bility of the numerical integration and the regularization method for infrared 
and collinear singularities, which needs to be performed in four dimensions. 

In the following sections we illustrate the method for the simple case of the 
(f) 3 model, showing some numerical comparisons between the results obtained 
with this approach and the standard diagrammatic techniques. 
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2 Introduction to the method 



We shall first shortly review the largest time equation and the dispersion rela- 
tions "a la Veltman" we shall use in this paper, referring to [25] for furthers^! 
details. 

For definitness we work in a 3 theory, setting the coupling constant to one 
for simplicity. The propagator for a scalar field, propagating from spacetime 
point x to y, is 

4* = / ^""p^ = *(*o - */o)A+ + % - *o)A+ (1) 

A ^ = / J^ e ~ U(X ~ y)9{1 ^ 12 ~ m ') ( 2 ) 

For a generic n-point, one-particle irreducible, Feynman diagram G n we obtain, 
in configuration space: 

Gn^Xi, . . . X n ) i A X \X2^X2%3 ■ ■ ■ Ax n a:i (3) 



where Xj enters clockwise in the diagram (see fig. d]). We now introduce a new 
notation: the "underlined configuration variable" x_j and a new set of functions 

Gn (%i i • • • i -ILj j • • • j 3Lki • • • ^n) 

(-lJTA^ . . . A. Xj _ lx _. . . . A XkX _ k+i . . . A,,.„ (4) 

where n u is the number of underlyined variables and 





= A + 




= A+ 




= A* 



Notice that because of the 9(Iq) term in equation (j2j) energy is alway required 
to flow towards the underlined vertex. 



1 We follow closely ref. [25], using Minkowski metrics and this justify some sign 
difference and the appearence of constant factors i. We also use the relationship 
A+ = A~ everywhere in the paper. 



4 




Fig. 1. Four point function Gi{x\, X2, £3, x^). 

We can introduce a diagrammatic notation for the G n functions: we shall de- 
note with a small circle a vertex to which an underlined variable is "attached" , 
namely, for example: 



x < 



G(x, y, z) = 



Therefore the set of Feynman rules for an underlined diagram are: 



x 



= A 



x 



•v = A+ • 



x 



= A* • 

xyi 



x 



A+ 

yx 



With the above set of Feynman rules and definitions we obtain 

x > y , z => G(x,y, z) = -G(x,y,z); G(x,y, z) = -G(x,y, z); . . . (5) 



namely if Xq is the largest time "entering into a given graph" the sum of a 
G with underlined x and without underlined x is equal to zero. This feature 
applies to any diagram with an arbitrary number of legs. This observation 
leads to the largest time equation 

G(x 1 ,...,x n ) = (6) 

all possible underlings 
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which, for the case of a three point function, reads 



G(x, y, z) + G(x, y, z) + G(x, y, z) + G(x, y, z) + G(x, y, z) + G(x, y, z) 
+G(x,y,z) + G(x,y,z) = (7) 

Indeed if x > y , z the contributions (1,2), (3,5), (4,6) and (7,8) of Eqn. (J7J) 
cancel pairwise. If y > x , z the cancellation occurs among the (1,3), (2,5), 
(4,7) and (6,8) pair of diagrams. Finally if z > x ,yo the cancellation occurs 
among the (1,4), (2,6), (3,7) and (4,8) pair of diagrams. This cancellation is 
explicitly worked out in the appendix. From Eqn. ([7]) one easily arrive at the 
standard Cutkowsky rules for the imaginarjd] part of a diagram 



G(x, y, z) + G(x, y, z) = G(x, y, z) + G(x, y, z)* 

= - G(x, y, z) + G(x, y, z) + G(x, y, z) + G(x, y, z) 
+G(x,y,z) + G(x,y,z) 



Since the Feynman rule for a line connecting an ordinary vertex and an un- 
derlined one is indeed equivalent to cut the diagram on the corresponding line 
this is the standard Cutkowsky rule for the imaginary part of the diagram. 
Notice that several contributions in the above equation vanish, since energy 
is required to flow from standard to underlined vertices and conflicting condi- 
tions might be required to hold simultaneously. By mean of the largest time 
equation we can derive another equation which allows to formulate a universal 
dispersive like formula for an arbitrary Feynman graph. We choose anyT^I two 
vertices in the diagram of interest, say y and z in the three point function. We 
have 

G(x, y, z) = 6(y - z )G(x, y, z) + 9(z - y )G(x, y, z) (9) 



and 



%o - z ) G(x, y, z) + G(x, y, z) + G(x, y, z) + G(x, y, z) 



(10) 



Indeed, as discussed in Eqn. ([5]), if y > z ,x contributions (1,2) and (3,4) 
cancel pairwise; if Xq > yo,z contributions (1,3) and (2,4) cancel pairwise; if 

2 Notice that, using the conventional decomposition of the S matrix S = I + iT, 
G(x±, . . . , x n ) contributes to iT 

3 that's why we choose to discuss both the simpler two point function and the 
three point one: for the two point function the arbitrairness in the choice of the two 
vertices has no meaning. 
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zq > Xo,yo the overall 9 function ensures the vanishing of this contribution. 
From Eqns. (|7j), Q and (ITU1) (and the analogous one for the term with 9(z —y ) 
we obtain 



G{x, y, z) = -6(y - z ) [G(x, y, z) + G(x, y, z) + G(x, y, z) 
-9(z - y ) \G(x, y, z) + G(x, y, z) + G(x, y, z) 
= -9(ya - zq) [G(x, y, z) + G(x, y, z) 
-O(zo-yo) G(x,y,z) + G(x,y,z) -G(x,y,z) 



Notice that the last term in Eqn. (TTTT) corresponds to a cut in which the two 
chosen vertices y and z lie on the same side of the cut. This term is one of the 
terms occurring in the unitarity relation (jSJ) and is purely real. 

For twice the imaginary part of G (real part of the amplitude), from Eqn. tjTTT) . 
we obtain 



G(x, y, z) - G(x, y, z)* = [6(y - z ) - 6(z - y )] G(x, y, z) + G(x, y, z) 



-G(x,y 7 z) -G(x,y, z) 
For the two point function we obtain 



(12) 



G(x, y) - G(x, y)* = [6(x - y ) - 0(y - x )] \G(x, y) + G(y, x) 
G(x, y) + G(x, yf = \G(x, y) + G{y, x 



(13) 
(14) 



This is a dispersion relation for the real part of a generic, one particle irre- 
ducible n-point function. In the appendix we verify explicitly that Eqn. (TT3T) 
holds and we show that inserting the integral representation of the 9 function 
it leads to a dispersive like integral. 

We introduce a diagrammatic representation for the 9 function appearing in 
Eqns. f|T2l) and f|T3|) : we will connect the two chosen vertices with a wavy line 
herafter referred to as the r-line. Before providing the generic recipe, notice 
that a function G is non vanishing only if ordinary and underlined vertices 
belong to two distinct connected regions: the one containing underlined vari- 
ables will be referred to as shadowed the other one as unshadowed. The energy 
along the cut lines will flow from the unshaded to the shaded region. We refer 
to [25] for a thorough discussion and we just give an example of a vanishing 
contribution to a box function in fig. [2j 

For a generic Green function the method can be applied as well and, sum- 
marizing the results sketched above, leads to the following modified Feynman 
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Pa Ps 

Fig. 2. Vanishing contribution to the four point function. The diagram 
is split into four disconnected regions (each propagator connecting a 
dotted and an undotted vertex is on shell) . Arrows denote the direction 
of the energy flow 

rules: 

• Draw a r— line (the 9 function insertion in Eqn. between two arbitrary 
vertices of the Feynman graph. 

• Cut the graph in two disconnected parts in such a way that the r— line 
crosses (if both the chosen vertices are on the same side of the cut this will 
contribute only to the imaginary part of the graph as G(x, y, z) in Eqn . (II II) ) 
the cut. 

• In the vertices where the r— line is absorbed/emitted it does contribute (this 
and the following items are esplicitly shown, for the two point function, in 
the appendix) to the energy balance with a r. 

• The r— line contributes a factor l/(«vr) x j dr (1/r) (f denotes the principal 
value integral). 

• Cut internal lines contribute a factor 6(±k°)5(k 2 — m 2 ) where ±k° is chosen 
in such a way that their energy flows towards the shaded region. 

• A two-body phase space integral is associated with the two cut internal lines: 
rf$ = \ l / 2 (P 2 ,p\,p%) / P 2 dcos9d(fi where pi2 are the momenta of the cut 
internal lines, and P is the sum of the external 4-momenta kj which enter 
into the unshaded part of the graph plus the r— line contribution which 
amounts to subtract r to the time-like P component. 

• Sum over the contributions of all allowed cuts (see the items above). 



2. 1 The two point function 

In order to illustrate the method let us start by computing the simplest 1-loop 
diagram in a scalar 3 theory (q is off-shell): 
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q-p 

Q 




P 

Using the standard techniques and omitting all constant factors for the sake 
of simplicity, this loop diagram is proportional to: 

/ dip {p*-m* + ie){{ q -py + ie) = / ^ / ^ V + L{x)Y ^ 



where one of the fields in the loop is assumed to be massless and 

L(x) = x(l — x)q 2 — (1 — x)m 2 . (16) 



Upon rotation to euclidean momenta and setting an ultraviolet cut-off A, the 
loop integral in (IT5l) reduces to 



rr 2 i dx dz- — — - ~ 

Jo Jo - L(x)r 

~ 7c 2 i I dx[ln A 2 - ln(-x(l - x)q 2 + (1 - x)m 2 )} = 
Jo 

2- fi A 2 (J 2 TTL 2 2 21 m% 1 2 , o 

= 11 InA In 5 -m — mm +2 

where in the latter step we have used integration by parts, and we have also 
used A 2 > q 2 ,m 2 . 

The largest time equation (|14p . for the two poin function, reads 



G(x,y)+G(x, y y = (A: y y + (A+ x y 

G(q) + G(q)* = -2^ q -^-6{q 2 - m 2 ) 

G(q) denoting the fourier transform of G(x,y). 
In fact, we have: 
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q-p 




where the parenthesis '[' indicates that the shadowed part of the diagram is 
on the right and vice- versa for ']'. The rule for internal lines flowing towards 
the shaded region is to associate to them a 8(k°)5(k 2 — m). For a line flowing 
out of the shaded region one associates a 8(—k°)5(k 2 — m). As a result, the 
latter diagrammatic equation reduces to (apart from an overall sign) 

j d 4 p[6(-q° + p°)6((q - p) 2 )6(p°)6(p 2 - m 2 ) + 

+6(q° - p°)6((q - p) 2 )9(-p°)8(p 2 - m 2 )} (17) 

The first term of the latter expression can be written introducing a Dirac-5 as 
j d A p J d A Q5 A (q-p-Q)6(Q )5(Q 2 )6(p°)5(p 2 -m 2 ). (18) 

which simply reduces to the phase space factor: 

f d 3 pd 3 Q 4 7i \ 1/2 (g 2 ,m 2 ,0) 7r q 2 -m 2 2 2 
J WW 5 {q ~ p - Q)= 2 q 2 = 2 -^ 9(q - m } ' (19) 

where the 9 factor is the threshold condition q 2 > m 2 . For simplicity, we have 
neglected constant factors. Computing also the second integral and accounting 
for the proper factors of n, one finds (see eq. ( fTTjl ): 

F + F* = -2n 3 q2 ~f 9(q 2 - m 2 ) (20) 



Let us now turn to the remaining part of our loop diagram, namely the one 
obtainable by subtracting F — F* . For this purpose, following |25j we consider 
the diagrammatic equation: 
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q-p 




Here we have a new symbol, the so called r-line represented by a / "5555W\ 
The r-line makes energy flow from the unshadowed to the shadowed areas Pi . 
Namely the 4-momentum associated to a r-line is r M = (r, 0,0,0). Following 
the prescriptions obtained in [25], the cut r-line introduces a factor 1/m x 
1/r and a principal value integral is required. Then the latter diagrammatic 
equation writes to (apart from overall constant factors): 



m J t [2 (q — r) 2 

7r (q + r) 2 — m 2 



6{{q - rf - m 2 ) 



1 f dr \ 7r (q — r) 2 — m 2 

L J _ b 

2 e((q + rf-m 2 )\. (21) 



2 (g + r) 5 

Let us make the shifts q°—r = —t' — ► r in the first integral and q°+r = r' — >• r 
in the second to obtain 



, 1 r 2 -q 2 -m 2 2 

T7 ^ 2 ^ r _c l + 

g u + r r z — q z 



1 r 2 — q 2 — m 2 
t — q° r 2 — q 2 



6(t 2 - q 2 - m 2 ) . (22) 



Apart from the 1/r factor, the cut diagrams are invariant and may be eval 
uated in any frame: we go to the frame q = 0, where the integral in (1221 



4 More precisely the T-line ensures: 1) that the energy is conserved on both sides 
of the cut and 2) that the cut lines exit from the unshaded area and enter into 
the shadowed one. Therefore the r-line accounts for an arbitrary amount of energy 
entering (most likely) the unshadowed area and exiting from the shadowed one, as 
indicated by the arrows in the figure. 
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reduces to 

* Ly^V -m')}. (23) 



Setting t' = —r and r = A we find 
A2 , A-m 2 1 

l2 dA — w- ( 24 ) 



Solving this integral and using the fact that A 2 ^> g 2 we find: 

F — F* = In A 2 — g2 " m ' ln(m 2 - q 2 ) - ^ lnm 2 . (25) 

g 2 g 2 



With the proper constants taken into account we reconstruct the result ob- 
tained initialled]- Indeed from our initial calculation we get: 



F - F* = 2n 2 i fin A 2 - ^ ln(m 2 - g 2 ) - ^ lnm 2 ) . (26) 

V Q ' Q J 




Fig. 3. Three point Green function in a scalar c/> 3 theory 
2.2 The three point function 



We consider now a three point scalar function in 3 theory. The Green function 
depicted in Fig. [3] is, up to constant factors, 



G 3 (P 1 ,P 2 ,P 3 )^ J d 4 L 



1 



[(L - A) 2 - m 2 ] [(L - Pi - P 2 ) 2 - m 2 ] [L 2 - m 2 ] ' 



5 Actually the two results differ by a constant factor. This is a consequence of the 
fact that the chosen Green function is divergent and the choice of two different UV 
regulators leads to a constant shift in the lagrangian parameters. 
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where Pj = (Uj,pj) denotes external momenta, all momenta are assumed 
incoming, m is the mass of the scalar particle and L = (L , 1) is the loop four 
momentum and the standard L — > L + ie prescription is understood. 

Using the rules explained above we can write: 



G 3 (Pi, P 2 , P 3 ) ~A 1 + A 2 + A 3 + A 4 , 

Ax~ [ f dk 1 dk 2 5(Tl 1 - r - E 1 - E 2 

J r + ie J 

5(pi-ki-k 2 ) 



E X E 2 \{K X + P 3 ) 2 - m 2 ] ' 

A 2 ~ / — T — [ dk 1 dk 2 5(n 2 + n 3 - r - E 1 - E 2 ) 
J T + ieJ 

S(p 2 + p 3 - ki - k 2 ) — , 

PiP 2 [(JCi - P 3 ) 2 - m 2 J 

^3 ~ / / dk 1 dk 2 5(n 1 + n 3 - r - E 1 - P 2 ) 

J r + ie J 

#(Pi +P3 - ki - k 2 )- 



P X P 2 [(^ - P 3 ) 2 - m 2 



A 4 ~ / / dkidk 2 <J(n 2 — t — E\— E 2 ) 



r + ie 



5(p 2 -ki-k 2 ) 



E X E 2 \{K X + P 3 ) 2 - m 2 



where Kj = (P,-, kj) = (y |kj| 2 + m 2 , kj) 

We can draw a r-line between Pi and P 2 vertices. We then cut (put on mass- 
shell) the corresponding internal lines and one of the other internal lines (this 
gives rise to two contributions corresponding to Ej both positive or negative). 
We sum over all possible cuts and convolute with the given weight as shown 
in Fig. HI 

Notice that we have obtained a dispersive formula: each term is the convolution 
of two tree level Feynman diagrams (with on-shell external particles and with a 
peculiar vertex where energy conservation is guaranteed only by including the 
r— line contribution) times two body phase space times 1/r dr integration very 
much like the dispersive integral of the imaginary part of an amplitude. This 
particular dispersive formulation is independent from the specific diagram. 

We can go one step further. Let us focus on the contribution A\\ it is the 
convolution of the process 
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Fig. 4. Three point Green function in a scalar cf) 3 theory. / "5MW v line 
denotes the dispersive integral over dr as well as internal two body 
phase space. The arrows on internal cut lines depict the energy flow in 
the cut diagram. 




Fig. 5. Convolution of amplitudes rather than individual Feynman 
graphs. It is manifest that, in addition to vertex contributions, we pick 
up also self energies contributions. 

and the t-channel graph for the process 

^cut) ^cut) ^ hh ^ (27) 
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where <pj denotes external particles with momentum Pj and (J)j denotes 
internal (cut) particles. If we extend A\ to include the contribution of the 
whole amplitude (s, t, and u channel contributions, see Fig. EJ) of process (127)) . 
we also obtain the contribution of the <\>\ external self-energies as well as part 
of the triangles for the process 

01 -> 0203 • (28) 



Applying the same procedure to all the A, graphs, we obtain the one-loop 
correction to the amplitude for the process (128)) as given by the formula 



0i - 2 03 = E (M» - 0^0^) * (<f>[ cut) 4 cut) - 0, 3 ) 

M>J2.J3 =1 » 3 

+ (0,, - * (0i CUt) 02 Cnt) - 0i 2 0i 3 ) , 

where with the symbol * we denote the convolution with the two body phase 
space and the dispersive integral over dr. Actually, to define unambiguously 
the above expression one needs a prescription to pick up the two vertices 
where energy conservation is ensured by the r-line contribution, and one needs 
to check that each diagram is accounted for with the correct combinatorial 
factors. This will be discussed later on. 

We have thus achieved our first goal: with the algorithm sketched above we 
can write any one loop amplitude as a convolution in 4-dimensions of tree level 
amplitudes. More precisely, given a set of external particles, we divide it into 
two non empty subsets Q a . j = 1,2. We compute the convolution of 

{G ai -> X1X2) * (X1X2 -> G a2 ) , (29) 



summing over all possible intermediate particles Xj (here "possible" means 
compatible with the symmetries of the theory). Recall that tree-level ampli- 
tudes are computed for on-shell particles but an arbitrary amount of energy 
can outflow/sink from/in one vertex. Therefore one of the two amplitudes in 
( 1291) can be non zero also for all particles incoming/outgoing. Finally we sum 
over all possible G ai , Q a2 partitions. 



3 The general algorithm 



In this section we outline the general algorithm we use to draw the r-lines 
and to select all the cut diagrams with the proper weights and combinatorial 
factors. 
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For definiteness we consider the amplitude 

0(Px)0(P 2 ) - 0(P 3 )0(P 4 ) (30) 

in (f) 3 theory. Pj denote, as before, the momenta of the external particles. 

• We divide the set of external particles into two non empty subsets Q aj 

j = 1,2. There are 14 f= Si=i w JJ sucn partitions. 

• We compute the convolution of 

[g ai - <p(K 1 )<p(K 2 )} * [0(^)0(^2) - g a2 ] (3i) 

summing over all possible partitions. We have denoted with Kj the momenta 
of the cut lines. 

• Require that the r line is emitted/absorbed at the opposite ends of the 
cut-line carrying K\ momentum. 

• Select an external momentum (arbitrarily). For definiteness we choose P\. 
We will call it reference momentum from here on. 

• If the reference momentum belongs to g a2 the amplitude for 
M. 2 = 4>{Ki)(j)(K2) — > g a2 needs to be modified as follows: 

• If a diagram has a propagator depending on a momentum Q = J2jQj, 
where Qj ^ K\ and at least Pi and K 2 enter the sum over j, this diagram 
must be vetoed. 

• If a diagram has a propagator containing K\ and at least an external Pj 
but not Pi, this diagram must be vetoed. 

■ The diagram of Ai 2 with K\ and K 2 attached to the same vertex should 
be vetoed if g ai is made up of a single external particle. This prescrip- 
tion allows us to avoid to compute external particle self energies. This 
contributions will be computed analytically and added in a second step. 

• The above prescriptions will apply also to g ai — > (p(Ki)(p(K 2 ) if Pi belongs 

to Gai ■ 

Let us examine one of such partitions: 

Gax = {0(^)1; g*2 = MPl), 0(A), 0(P 3 )} 

The calculation of the convolution M. ai * M. a2 f° r such partition is sketched 
in Fig. [TDJ Vetoed Feynman diagrams are marked in the figure. Summing up 
over all fourteen partitions one can check by inspection that all contributions 
to NLO correction are correctly taken into account and that any relevant cut 
occurs the correct number of times. 

After providing the recipe, let us try to motivate it. We shall discuss one of the 
box diagrams contributing to the amplitude in Eqn. (150]) . As shown in fig. E] 
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there are six possible different ways to draw the r-line, four corresponding to 
the r-line emitted and absorbed by two nearby vertices and two corresponding 
to the r-line emitted and absorbed by two opposite vertices. 



Pi 

--'C4 






. Pi 


Pi 






P 2 






% 






Cs 






Pa 




p 3 


Pa 




^3 



Fig. 6. Four point function. There are six possible ways to draw the 
r-line: four as shown in the left-hand side of the figure (the other three 
possibilities are similar but with the r-line among P2 and P3 or P3 and 
P4 or P4 and Pi) and two as shown in the right-hand side of the figure 
(the other one is similar but with the r-line among P2 and P4). The left 
hand choice gives the sum of six contributions corresponding to having 
the propagator connecting Pi and P2 cut together with anyone of the 
remaining propagators. The right hand choice gives eight contributions 
corresponding to having a cut propagator between Pi and P2 or be- 
tween P2 and P3, together with the propagator between P3 and P4 or 
between P4 and Pi. Recall that each "cut" actually corresponds to two 
contributions depending on the energy flow. 

Let us focus on the case of the r-line between two near-by vertices: the 
diagrams will be written as the sum over six possible cuts. If we denote by 
TTij the propagator connecting Pj and Pj and by (71^, 7T/ m ) the contribution 
corresponding to iij k and 7i [rn cut propagators, the left-hand diagram of fig. [6] 
will result from the sum of the (hu^ai), (tti2, ^34) and (7Ti2,7r 2 3) cuts, each 
counting as two distinct cuts dependig on the energy flow. Each of these con- 
tributions will be included into a different term of the convolution of Eqn. (13ip . 
namely it will arise with a different choice of otj (notice that on the opposite 
sides of the cuts in Fig. [6] there are different sets of external particles). There- 
fore we need a prescription to ensure that all the cuts enter into the final sum 
with the same combinatorial factor and that we do not have contributions 
from the other possible ways of drawing the r-line. 

The simplest prescription is to require that the r-line is emitted and absorbed 
together with particle 4>(Ki), namely it will enter/exit always from the vertices 



6 The other choice is possible as well. It requires a somewhat more involved algo- 
rithm. We have not studied it insofar since at present we do not have a good reason 
to believe it leads to easier integrals. 
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to which 4>(Ki) is attached. In this way the r-line is always emitted and reab- 
sorbed at the opposite side of one propagator and we single out the topology 
we are after (left panel of Fig. [HD- 

If we allow 4>{Ki) to connect to all possible external particles, the diagram of 
Fig. [6] (left panel) would be reconstructed four times, since the r-line would 
connect all possible near-by vertices. With the same reasoning it's easy to 
conclude that a n-point diagram will be recovered n times. 

To adjust the combinatorial factors, we introduce the reference momentum 
Pi: we require that <f){Ki) is attached to the same vertex with <fi(Pi) or with 
a tree of external particles containing 0(Pi) as shown in fig. [7J Notice that 
in this way we obtain each diagram twice (the <fi{Ki) vertex has two near-by 
vertices) with the exception of self-energy contributions which are generated 
only once. These are indeed the correct combinatorial factors. 



Pr 



Pi 



K x P2 



P 3 K x P 2 



Fig. 7. Allowed particles tree into the amplitudes entering into the con- 
volution of Eqn. (f3T|) . Cut propagators (denoted with a small circle) 
carry momentum Kj and external particles carry momentum Pj. The 
r-line is always emitted/absorbed from a vertex connected either di- 
rectly to Pi (left-hand side) or to a tree made up of external particles 
only (K2 not allowed) and containing Pi . Equivalently, if the tree is seen 
from the opposite side, the r-line is always emitted/absorbed from a 
vertex connected to a tree made up of external particles plus K2 and 
not containing Pi. The left-hand side diagram, inside the convolution, 
will contribute to the evaluation of boxes and triangles, whereas the 
right-hand side to triangles only. 

All the prescriptions given above to veto certain contributions to the ampli- 
tudes aim at implementing the above reasoning. They amount to require that 
any particle tree starting from 4>(Pi), to be accepted, should absorbe 4>(Ki) be- 
fore 0(^2), as shown in Figs. (J7J) and (jBJ). They look somewhat involved since, 
due to four-momentum conservation, the momentum flowing into a given prop- 
agator can be computed via two possible combinations of particles and we do 
not know a priori which will be chosen from our tool for the automatic evalu- 
ation of tree-level amplitudes. Therefore we cannot simply check if K\ and Pi 
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enter into the "same" (in the extended sense described above) vertex, but we 
need also to check (and accept) also trees made up of external particles 0(Pj) 
{j ^ 1) and/or <f>{K 2 ) 



K x P x 



Ki P 3 P 1 



Fig. 8. Vetoed particles tree into the amplitudes entering into the con- 
volution of Eqn. (I3ip . Cut propagators (denoted with a small circle) 
carry momentum Kj and external particles carry momentum Pj. The 
K 2 particle is always emitted/absorbed from a vertex connected either 
directly to Pi (left-hand side) or to a tree made up of external parti- 
cle only (if 1 not allowed) and containing Pi . Equivalently, if the tree is 
seen from the opposite side, the K x particle is always emitted/absorbed 
from a vertex connected to a tree, made up of external particles, which 
does not contain neither K 2 nor Pi. 

Coming back to the four point function we sketch in fig. [9] the way it arises 
from the convolution in Eqn. (I3TI) . For simplicity we limit to a single box 
contributing to the amplitude and denote with dots all the remaining contri- 
butions. As shown in fig. [9] the box is evaluated twice corresponding to two 
different ways of drawing the r-line. 

The various cuts, denoted with Cj in fig. [9j arise as follows in Eqn. (I3T1) 



[01- 


(j>{K x )(j>{K 2 )] * 


[<f>{K x )<f>{K 2 ) - 


02,3,4] 


=c 2 + c 8 + 


[02- 


<t>{K x )<t>{K 2 )\ * 


[4>(K X )4>(K 2 ) - 


01,3,4] 


=c 6 + ... 


[0 3 - 


(j>{K x )(j>{K 2 )\ * 


14>(K X )4>(K 2 ) - 


01,2,4] 




[04- 


(j>{K x )(j){K 2 )\ * 


14>(K X )4>(K 2 ) - 


01,2,3] 


=c 12 + ... 


[^1,2 - 


- 0(^i)0(^ 2 )] 


* [0(^)0(^2) - 


-» 03,4] 


= Cio + ... 


[01,3 - 


- 0(^i)0(^ 2 )] 


* [0(^)0(^2) - 


02,4] 




[01,4 - 


-> 0(^)0(^2)] 


* [0(^)0(^2) - 


-* 02,3] 


=c 4 + ... 


[02,3 - 


- <f>(K X )cf>(K 2 )} 


* [0(^)0(^2) - 


- 01,4] 


=c 3 + ... 


[02,4 - 


-> <f>(K X )cf>(K 2 )} 


* [0(^)0(^2) - 


- 01,3] 




[03,4 - 


- 0(^l)0(^2)] 


* [0(^)0(^2) - 


- 01,2] 


=c 9 + ... 


[01,2,3 


- (f>{K x )cf>{K 2 ) 


] * [<j>{K x )(j>{K 2 ) 


-04] 


=c lx + ... 


[01,2,4 


- 4>(K X )4>(K 2 ) 


] * [<P{K X )<P{K 2 ) 


-03] 




[01,3,4 


- (f>{K x )cf>{K 2 ) 


] * [<P{K X )<P{K 2 ) 


-02] 


=c 5 + ... 


[02,3,4 


- (f>(K x )(f>{K 2 ) 


] * [<P{K X )<P{K 2 ) 


-0l] 


=d + c 7 .. 
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= £,•[&„ - 0(ATi)0(/f 2 )] * [*(A'i)*(A- 2 ) -> ft,,] 



4 




C7/ Cg 




+ 



Fig. 9. Cuts contributing to the four point function. Only those relevant 
to a specific box are shown 



The above algorithm applies also to more complex cases with one important 
addition: 



• If cut particles are not identical, one needs to consider both the contribu- 
tion with K\ <-> K2. Namely the contribution of r-line absorbed/emitted 
from K\ with the above vetoing prescriptions plus the contribution of r-line 
absorbed/emitted from K 2 with the above vetoing prescriptions but with 
Ki <-> K 2 . The sum of this two contributions leads to twice the results 
and again one can check by direct inspection that the combinatorics of cut 
diagrams is correct. 

The entire procedure is numerical: we perform a numerical integration si- 
multaneously over the phase space variables of external particles, over dr and 
phase space variables of intermediate particles. Each computation is performed 
in four space-time dimensions. If needed the integral in dr is truncated at 
r < T max (where r max is an ultraviolet cutoff) and proper counterterms are in- 
cluded to ensure that the integral is properly regularized. In practice one com- 
putes a subtracted Green function and then modifies the tree-level Lagrangian 
with appropriate counterterms in such a way to reproduce the correct results 
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for divergent Green functions. 




-Pi P3 P2 Pi P3 P2 

Fig. 10. Convolution of amplitudes rather than individual Feynman 
graphs. Cut lines carry an arrow. / WSSW S indicates a r-line. The graphs 
contained in the picture correspond to one of the fourteen partitions of 
the external momenta. The non-vetoed graphs correspond to different 
one- loop diagrams: ai, £14 and 05 are cuts contributing to three different 
triangles, whereas 02 and contribute to two different box diagrams. 
Considering the whole set of fourteen partitions, one recovers the full 
four-point Green function. 
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4 Numerical results 



We give, in this section, few examples of results obtained with the above 
described numerical method compared to standard analytic evaluation of loop 
integrals More complete numerical results will be shown elsewhere. First 
of all we give the predictions in Tabs. [T], [2J [3] for the three-, four- and five- 
points scalar form factors respectively, defined in Eq. [321 in the <f) 3 scalar 
theory. We compare only the real part of the form factors, since the method 
illustrated in the previous sections allows to calculate it. The imaginary part 
can be easily calculated with standard Cutkowsky rules. This is much easier 
than the real part, since the required integration over the internal cut lines is 
only the 2-body phase-space and no dispersion integral is required. Moreover 
no principal value integration is neede 



If 1 

Co = — / d q , 

D = —fd A q- 1 



.spei 



in 2 J 6/1^2^3(^4' 

E = ^ [ d 4 q ) , (32) 
m z J d^d^d^dr, 



where 



di = (q 2 — m 2 + ie) 

7 For the real part of the amplitude we deal with the numerical evaluation of 
principal value integrals using the following strategy. For simplicity let's assume we 
have to perform the principal value integral of a function f(x) with a pole in xo- 
We can either perform the integral of 

/!(x) = i[/(x) + /(2x -x)] 
or, use integration by part, 

/ f(x)- — —dx -> log|x - x Q \(x - x )f(x) 

J X-Xq 

/d 
log \x - x o\-^ii x - x o)f( x )]dx 

in both cases obtaining a numerically convergent integral. Both methods have been 
explored and tested. (We however need to study computationally more involved 
examples to be sure that numerical instabilities do not arise because of nearby 
poles or soft/collinear enhancements). Notice that we can apply the method here 
sketched since we can easily evaluate analytically the location of the poles as a 
function of the integration variables. 



22 





^ ( t\ 
p x (CxeVJ 


Py (GeV) 


p 2 (GeV) 


Pi 


1. 1 ozUOUol 


1 /I1/I91Q^/I 
1.414Z1O04 


— Z.ZoDUDoUl 


P2 


-1.73205081 


-1.41421354 


2.23606801 


P3 













Re(C ) 


1.982(2) (numerical) 


1.98390995 Ref. [29] 



Table 1 

Comparison between the values of the scalar form factor Co for the <^> 3 theory as ob- 
tained with the numerical method and with the analytical results of LoopTools |29j . 
The momenta components are taken to be p = (Px,Py,Pz)i with a common internal 
and external mass of 0.01 GeV. 





p x (GeV) 


Py (GeV) 


Pz (GeV) 


Pi 


1.73205081 


1.41421354 


-2.23606801 


P2 


-1.73205081 


-1.41421354 


2.23606801 


P3 


-2.23606801 


-1.41421354 


-1.73205081 


P4 


2.23606801 


1.41421354 


1.73205081 




Re(A)) 


0.4435(2) (numerical) 


0.443615191 Ref. [29] 



Table 2 

Comparison between the values of the scalar form factor Dq for the <p 3 theory as ob- 
tained with the numerical method and with the analytical results of LoopTools [29] . 
The momenta components are taken to be p = (Px,Py,Pz), with a common internal 
and external mass of 0.01 GeV. 





Px (GeV) 


Py (GeV) 


Pz (GeV) 


Pi 








2.13816766718398 


P2 








-2.13816766718398 


P3 


-0.825139760971069 


-0.878521561622620 


-0.08667117357254028 


P4 


-0.501121819019318 


-0.772821187973022 


-0.162467777729034 


PS 


1.32626157999039 


1.65134274959564 


0.249138951301575 




Re(£ ) 


-10.48(2) (numerical) 


-10.4724461 Ref. [29] 



Table 3 

Comparison between the values of the scalar form factors for the cp 3 theory as ob- 
tained with the numerical method and with the analytical results of LoopTools [29] . 
The momenta components are taken to be p = (p x ,Py,Pz), with a common internal 
and external mass of 0.01 GeV. 

di = ((<? + h-i) 2 -m 2 + ie) 
and hi are related to the external momenta as pi — k\, p n — k n — fc( n -i)- 
In addition we compare the pure virtual NLO QCD correction to e + e _ — > 
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1611111011 llldbb V J 


dlldiyilCdl ^IIU J 


nmiieilCdl I ilU J 


1 
1 


-l.U0ZZ(^Z J 


-1 .UOZO^O J 


Z 


-U. lo l\)[L ) 


-U. /ooz^o ) 


3 


-0.5985(1) 


-0.5989(2) 


8 


+0.01320(0) 


+0.01289(8) 


9 


+0.10271(0) 


+0.10258(6) 



Table 4 

Values of the virtual QCD cross section for the process e + e~ — ► 7* — ► gg, at y^i = 
20 GeV, for different fermion masses. The infrared divergence has been regulated 
with a photon mass m 7 = 30 MeV. 

qq production (limited to the photon exchange case), at y/s = 20 GeV and 
for different values of the fermion mass. The infrared singularity has been 
regulated giving a mass of 30 MeV to the gluon. Since no QCD non-abelian 
vertex is present the gluon mass can be safely used. The QCD coupling g s has 
been set to 1 for simplicity. 

For e + e~ — > qq there is a divergent 3-point function. We deal with the renor- 
malization problem as follows: 

• We compute the 3-point function using a hard cut-off r max . 

• For a specific external momenta configuration we compute the 3-point Green 
function both using the hard cut-off r max , obtaining G^" 13 "^, and with a 
conventional regularization, e.g. dimensional regularization (DR), obtaining 

G (DR) 

• We then add to the Lagrangian the counterterm 

G (DR) _ G ( W ) = ^ (a (^ x ) 7oAo _ a ( Tmax ) 7 . 



5 Conclusions 

We have proposed a new approach to the evaluation of one-loop scattering 
amplitudes in perturbative quantum field theory. In the present paper we 
have presented the method via some basic examples of one-loop amplitudes in 
4>\ field theory. The basic rules of our calculation approach have been briefly 
illustrated. 

We have calculated few scalar Green functions and evaluated the QCD next- 
to-leading order correction to e + e~ — > 7* — > qq finding very good agreement 
with analytical results. 

To assess the viability of the approach there are still several steps to be faced. 
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We plan to: 



• Apply the method to the evaluation of n-photon amplitudes in QED. 

• Explore the feasibility of a regularization based on a hard momentum cut- 
off both in the ultraviolet and in the infrared region for non abelian gauge 
theories. This should be possible provided a suitable number of counterterms 
is introduced in the Lagrangian [27]. 

• Reproduce available results in e + e~ — ►jets at next-to-leading order accuracy. 

• Establish the relationship among our regularization scheme and the stan- 
dard dimensional regularization in order to be able to use existent parton 
distribution functions. 
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6 Appendix A 



A+- G ► z° > x°j 



\j • © > %i 



A* 3 G G 



Fig. 11. Feynman rules for propagators involving at least one dotted 
vertex 



For the self energy diagram we have (see Eqn.fJTJ)} 



G(x, y) = A xy A yx = (^A+ + )(t? wx A+ + $ xy A yx ) 



'xy ' " yx—-yxj \" yx^yx ' " xy^-yx 
■#xy(At y ) 2 + #yx(&ix) 2 

= 1? 



x$(k°)$(q°)5(k 2 - ml)6(q 2 - m 2 2 )} 
+ ^{(2tt)6 / d4fcd4g e w[~i(y ~ x)(k + q)] 
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X 



#(k°)#(q°)5(k 2 - m\)5{q 2 - m 2 2 )} , (33) 



where d yx = 9{xq — yo). We can now introduce an integral representation for 
the 9 function 

6{x) = — -^-dr, 



2m / t — it 



with e — > + and Eqn. (1331) becomes 

1 f exp[iT(x° — y )} 



r — xe 



2m J 

x |(2^)6 / d4/cd4 9 exp[-i(a;-2/)(A: + g)] 

x^(fc°)^(g°)5(A; 2 - m 2 )c% 2 - m 2 )} 



2m' J t — ie 



x i / d 3 A;d 3 g exp[— i(y — x)(A; + g)] 



x#(k )#(q°)6(k 2 - mj)5(q 2 - m 2 )} 



Since we want to compute the self energy in momentum space, we perform 
the Fourier transform 

G{p 1 ,p 2 ) = J e ipix e ip2y G(x,y) dx dy w= = +v J e ipix e ip2{w+x) G{w) dx dw 



and we obtain 



dr f d 4 k d 4 q—^—5 3 (p 2 + k + q)$(k )$(q )5(k 2 - m 2 )5(q 2 
J t — ie 

x[5(t + p 20 + k + q ) + 5(t - p 20 + k + q )) 



(2tt) 4 5(pi + p 2 ) / dr / d A k 



(2n) 3 i J r — ie J 2/c 2q 

x [5(t + p 20 + k + q ) + 5(t - p 20 + k + g )] 

where, to perform x and w integrations, we have used J e lkx d A x = (27r) 4 5 4 (/c) 
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and k = J k 2 + m 2 , q Q = J (p 2 + k) 2 + m 2 . Recalling that 



/(r) dr = / t^l dr + iTT / /(r)5(r) dr 



T — It 



we obtain 



G(pi,P2) = 77r ^r(27r) 4 5(p 1 +p 2 ) 




(27r) 3 r ' ^ \J r J 2k 2q 

x [S(t + p 2 o + h + g ) + S(t - P20 + ^0 + ?o)] 

+Z7T / d 3 /c-^--^(5(p 20 -ko- g ) [ (35) 
7 2/c 2g J 

where we assume p 20 > and thus S(p 20 + k + g ) = 0. We can work in the 



p 2 rest-frame where k = \jk 2 + m 2 , g = \Jk 2 + m 2 . The phase space volume 
d 3 A; can be written as k 2 dk dQ and the integral over dk can be performed 
eliminating the 5 functions in Eqn. (1351) and we finally obtain 



1 f p- m l- m 2-P2(l 1 f u 

G(p 1 ,p 2 ) = —X27ry6(p 1 +p 2 )U dr- j dQ 



-r - p 20 



+ + dr- dQ +Z7T / dQ^- 

r J ~r + p 20 J p 20 



k2 = (p 20 ± t) 2 _ m\ + m\ (mj - m 2 2 f 
± 4 2 4(p 20 ±r) 2 

fc2 = Po _ m 2 + m 2 + Oi - m 2 ) 2 



4 2 4^ 



The last term correspond to the real part of the diagram and is indeed equal to 
half the contribution found, in Eqn. fTSOj) . using the largest time equation. The 
first two terms correspond to the imaginary part of the integrals and coincide 
with the result of Eqn. (|25|) . Notice that the dispersive integral over dr arises 
as a consequence of using the integral repesentation for the 9 function, thus 
in the equation (fl3l) 



1 f 1 

[0(xo - yo) - 0(2/0 - x )\ -► — -/-- 

tlT J T 



and it is responsible for the appearence of the dispersive integral as well as of 
the appearence of r into energy momentum conservation 5 functions. 
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7 Appendix B 



We shall now explicitly verify Eqn. (j7]). We have 



(l/t 3 )G(x, y, z) = i 3 A xy A yz A zx = 9 xy 9 yz At y A+ A+ + 9 xy 9 zx A+ At A 



xy yz xz 1 " xy" zx" X y zy zx 

+9 zy 9 xz A x ~ y A+ y A+ z + 9 yx 9 xz A y ~ x A y ~ z A+ z 

i_n Q A + A + A + i $ $ A + A+ A + 
"r (7 j/2 (7 2;x^j / x' iA yz'- A z2: ~r u yx° zy L -±yx L -^ zy^zx 

9 9 A + A + A + 4- 9 9 A + A + A + 

u xy^yz'-^ X y'-^yz'-^xz \ y xy u ZX<-± x y<-± zy <-± zx 

+ A + A + A + - # A + A + A + 

^^xy^zy^xz v yz L -^xy L -^zy L ^-xz 

-9 ZT A + A+ A+ + 9 1IZ 9 ZX A + A+ A+ 

" zx^xy zy xz 1 ^ zy xz 

4- A+ A + A + - (9 A + A + A + 

^^j/x^yz^xz w xy L -*y x '- y yz L -*xz 

n A + A + A + _|_ fl A + A + A + 

C'zx^yx^i/z'^xz T V X yU zx i-iy X i-±y Z i-± xz 

+9 9 A + A + A + + A + A + A + —9 A + A + A + 

~ v yzV zx L -±y X L -±yz L -*-zx ' L -^yx L -^ zy 1 -* zx v xy L -^y X L -^zy L -^zx 

n a + a + a + i a a A+ A + A + 

u yz L -^y X L ^zy L -^zx ' u xy u yz L -^yx L ^zy L ^zx 

:A + A + A + + A+A+A+ + A+A+A+ 

xy zy xz 1 yx yz xz 1 yx zy zx 



-®yz [A X yA Z yA xz + \jx^zy^zx 

— 9 zx ( A + A+ A+ + A+ A+ A+ ) 

zx y—*xy zy xz yx yz xz J 

-9 ( A + A + A + 4- A + A + A + 1 

u xy \ L -*yx L -*yz l -*xz ' L -*yx L -*zy L -*zx) 

+9 xy 9 yz (A+ A+ A+ + A+A+ A 
+9 xy 9 zx (A+ A+ A+ + A+ A+ A 
+9 yz 9zx (A+A+ A+ + A+ A+ A 



where = #(^4o — -Bo) an d we have repeatedly used 6Ub = 1 — 
Analogously we have 



{l/i 3 )G(x,y,z 
(l/i 3 )G(x,y,z 

(l/i 3 )G(x,y,z 
(l/i^x^z 
{l/i 3 )G{x,y,z 
(l/i 3 )G(x,y,z 
(l/i 3 )G(x,y,z 



■ 9 11Z [ A + A + A + 

V z \ yx zx zy 

■ 9 Z x (A^ A+ A+ 



yx 

+ 
j/z 



^ (a+a+a 

■ 9 yz (A+A+A 
^ (A+A+A+ 
*«» (A+A+A+ 



A + A + A + 

yx zx yz 

A + A + A + 

XJ/ ZJ/ zx 



A + A + A + 

xz yz xy 



A+ A + A + 

yx zx zy 

A + A+ A+ 

xy zy xz 



A + A + A + 

xz yz yx 



A + A + A + ) + A + A + A + 

xy xz zy J 1 xy xz zy 



A + A + A + 

yx zy zx 

A + A + A + 

xz yz yx 



+ A + A + A + 

1 yx zy zx 

+ A + A + A + 

1 xz yz yx 



A+A+A+ -A+A+A+ -A+A+A+ 

yx yz zx xy zy zx xy yz xz 
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+6 V . ( A+ A+ A + + A + A+ A+ ) 

i "yz y—*yx yz zx ' xy yz xz J 

+6 zr ( A + A+ A + + A + A+ A+ ) 

i "zx y—*yx yz zx xy zy zx J 

+6 rv ( A+ A+ A+ + A+ A+ A+ ) 

i "xy \~xy zy zx ' xy yz xz J 

O n ( A + A + A + I A + A + A + A 

Vxy u yz \ L -^xy L -^yz L -^xz < LJ, -yx L -^zy L -^zx j 

-B„.B XT ( A + A + A + + A + A + A + ) 

w xy w zx y—^xyzyzx 1 yx yz xz J 

-6 yz e zx (A+A+ A+ + A+ A+ A+ ) (36) 
and from Eqns. 03611361) we indeed obtain Eqn.(J7J). 
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